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ABSTRACT 

Aims. The aim of this paper is to test the basic model of negative AGN feedback. According to this model, once the central black 
hole accretes at the Eddington limit and reaches a certain critical mass, AGN driven outflows blow out gas, suppressing star formation 
in the host galaxy and self-regulating black hole growth. 

Methods. We consider a sample of 224 quasars selected from the Sloan Digital Sky Survey (SDSS) at z < 1 observed in the infrared 
band by the Herschel Space Observatory in point source photometry mode. We evaluate the star formation rate in relation to several 
outflow signatures traced hy the [O III] 44959,5007 and [OII] 43726,3729 emission lines in about half of the sample with high quality 
spectra. 

Results. Most of the quasars show asymmetric and broad wings in [O III], which we interpret as outflow signatures. We separate the 
quasars in two groups, “weakly” and “strongly” outflowing, using three different criteria. When we compare the mean star formation 
rate in five redshift bins in the two groups, we find that the SFRs are comparable or slightly larger in the strongly outflowing quasars. 
We estimate the stellar mass from SED fitting and the quasars are distributed along the star formation main sequence, although with 
a large scatter. The scatter from this relation is uncorrelated with respect to the kinematic properties of the outflow. Moreover, for 
quasars dominated in the infrared by starburst or by AGN emission, we do not find any correlation between the star formation rate 
and the velocity of the outflow, a trend previously reported in the literature for pure starburst galaxies. 

Conclusions. We conclude that the basic AGN negative feedback scenario seems not to agree with our results. Although we use a 
large sample of quasars, we did not find any evidence that the star formation rate is suppressed in the presence of AGN driven outflows 
on large scale. A possibility is that feedback is effective over much longer timescales than those of single episodes of quasar activity. 

Key words. Galaxies:quasars; Galaxies: active; Galaxies: star formation. 


1. Introduction 

An important discovery in the last 15 years is that the black hole 
mass at the centre of the galaxies is a constant fraction of the 
mass of the stellar bulge (Ferrarese & Merritt 2000, Gebhardt 
et al. 2000). From that moment it has been clear that the evo¬ 
lution of galaxies is closely connected with the growth of the 
central black hole and that the influence of the black hole on the 
host galaxy cannot be ignored. 

In particular AGN feedback (i.e. the energy released by the 
AGN into the interstellar medium) is invoked to solve many as- 
trophysical problems. For example, the very sharp cut-off seen 
at the bright end of the galaxies luminosity function cannot be 
reproduced without invoking a feedback mechanism. In other 
words, super-winds may be responsible for the lack of bright 
galaxies that are not observed in reality (Benson et al. 2003) and 
for the red, dead elliptical galaxies we observe today. Another 
crucial problem that requires AGN feedback is the so-called 
cooling flow problem at the centre of galaxy clusters, which 
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requires a source of heating able to balance or quench cooling 
(Fabian 2012 and references therein). 

A wide variety of quenching mechanisms have been pro¬ 
posed to self-regulate the growth of the galaxy (i.e. the star for¬ 
mation rate) and the black hole mass accretion rate (e.g. Di Mat- 
teo et al. 2005; Hopkins et al. (2006); Croton 2006; Martig et al. 
2009). Generally, the quenching is identified with a triggering 
event, for example a major merger, that channels a large amount 
of gas into the central region in a short period of time, activat¬ 
ing both star formation (SF) and quasar activity. The resulting 
AGN feedback rapidly exhausts or sweeps away the gas from the 
galaxy. Theoretical models (e.g. Hopkins et al. 2006, Granato 
et al. 2004, Di Matteo et al. 2005) have proposed that feedback 
from AGNs is responsible for the regulation of SF, transforming 
a starburst galaxy into a red elliptical galaxy. 

In a considerable fraction of AGNs there are clear signa¬ 
tures of outflows likely able to remove significant amounts of 
cold gas from the galaxy. AGN-driven outflows with veloci¬ 
ties of ~ 1000 km s ' have been identified in ultra luminous in¬ 
fra red galaxies (ULIRGs), e.g. Cicone et al. 2013, Rupke & 
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Veilleux 2013, and fast outflows are seen in broad absorption 
line QSOs displaying broad UV absorption lines with widths of 
several thousand kms '(Turnshek et al. 1988). The highly ion¬ 
ized Fe absorption lines blueshifted by 0.05-0.3c observed in the 
X-ray bands (e.g. Tombesi et al. 2010) have been found asso¬ 
ciated with the much slower (~1000 km/s) outflow components 
seen in the molecular gas (e.g. Feruglio et al. 2015, Tombesi 
et al. 2015). On the galactic scale, integral held spectroscopic 
observations have revealed the presence of ionized gas outflows 
with velocites up to 1000 kms 'in low and high redshift galax¬ 
ies (Alexander et al. 2010, Harrison et al. 2012, Harrison et al. 
2014, Cano-Diaz et al. 2012). 

However, studies exploring the role of AGNs with respect 
to star formation have led to ambiguous and contrasting results, 
with evidence for both negative and positive impact of AGN out¬ 
flows on star formation. After many years of extensive research, 
conclusive proof that AGN feedback is able to halt galaxy-wide 
star formation is still lacking. Moreover, many studies show that 
the opposite may actually be true. The squeezing and compres¬ 
sion of cold gas induced by accretion related winds could cause 
local density fluctuations leading to star formation (Ishibashi & 
Fabian 2012). Theoretical models able to explain the connection 
between the central engine and the host galaxy properties with 
secular evolution without invoking AGN feedback have been de¬ 
veloped (e.g. Jahnke & Maccio 2011, Ciotti & Ostriker 2007). 
Finally, even if it is well known that the AGN driven outflows 
are common in both radio loud and radio quiet quasars, it is not 
obvious if there is an efficient coupling mechanism between the 
AGN and the ISM (e.g. Bicknell et al. 2000, Kalfountzou et al. 
2012 ). 

In this paper we focus on ionized gas outflows traced by the 
optical [OIII] and [Oil] lines. Since the forbidden line emis¬ 
sion is isotropic and self-absorption in narrow lines is negligi¬ 
ble, the blue wings are generally interpreted as the result of out¬ 
flowing gas. The [OIII] line is typical for AGNs and originates 
from the ionized narrow-line region (NLR) gas surrounding the 
accreting super massive black hole in the centre of the galaxy 
(see Osterbrock 1989). Mullaney et al. 2013 constructed aver¬ 
age [OIII]/15007 profiles derived from stacking SDSS spectra 
of type 1 AGNs and found profiles displaying prominent blue 
wings shifted with respect to the narrow components. Similar re¬ 
sults are found by Zakamska & Greene (2014) who analysed the 
SDSS spectra of obscured luminous quasars. The [OIII] 45007 
emission line typically shows blueshifts and blue excess, likely 
signatures that the NLR is undergoing an outflow. 

The question we address in this paper is whether SMBH out¬ 
flows and feedback have the ability to regulate the starburst on 
a galaxy-wide scale as predicted by the negative feedback sce¬ 
nario. To this end, we compare the kinematic of the ionized gas 
outflows derived from the [O III] line profile to the star formation 
rate of the host galaxies estimated through the fit of photo-metric 
points in the infrared. Since the two processes likely act on dif¬ 
ferent timescales (e.g. Stanley et al. 2015 and references therein) 
we need a large sample of luminous powerful quasars in order to 
recover the signature that AGN outflows affect the star forma¬ 
tion rate of the host galaxy. Direct Herschel FIR measurements 
of individual galaxies is a unique tool used to properly estimate 
the ongoing star formation rate, especially in low luminosity IR 
galaxies, that in large sky surveys would not be detected. There¬ 
fore, we selected a large sample of quasars at redshift < 1 (to 
cover the [O III] emission line wavelength in SDSS optical spec¬ 
tra) observed by the infrared telescope Herschel in photometric 
point source mode. 



Fig. 1. Histogram of the redshift distribution for the quasars in our 
sample. The red, blue, and black hatched areas identify quasars from 
the DR7Q, DR9Q, and SDSS archives, respectively. The yellow area 
denotes the sum of the three subsamples. 



Fig. 2. Herschel PACS and SPIRE images in different bands for a 
quasar from our sample. The green circle is centred on the source and 
represents the extraction region. 

The paper is organized as follows. In Sect. 2 we define the 
sample of quasars, and in Sect. 3 we present the analysis of the 
SDSS spectra and of the Hershel data. The main results and the 
comparison between the outflow properties, the SFR and AGN 
properties, are presented in Sects. 4 and 5. In Sect. 6 we provide 
a discussion, and in Sect. 7 our summary and conclusions. 

In this paper we assume the following cosmology parame¬ 
ters: Q.„-0.3, A = 0.7 and Ho = 70 km s“* Mpc“*. 

2. The sample 

The AGN parent sample includes bona fide quasars extracted 
from the Sloan Digital Sky Survey (SDSS) archive or the SDSS 
quasar catalogue that have high resolution photometric images 
from the infrared satellite Herschel. 

First of all, we selected a list of reliable quasars candidates 
from 

- the fifth edition of the spectroscopic quasar catalogue of 
Schneider et al. (2010) based on the SDSS seventh data re¬ 
lease (DR7Q). This catalogue contains 19552 quasars atred- 
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shifts of less than 0.8 (this redshift limit ensures that [O III] is 
in the observed frame). 

- the quasar catalogue of Paris et al. (2012) based on the SDSS 
ninth data release (DR9Q). There are 23964 quasars at red- 
shifts of less than 1. 

- the SDSS archive. We searched for objects classified as 
‘QSO’ in the latest SDSS data release (DRIO) not present in 
the two quasar catalogues with redshifts of less than 1 and a 
median signal-to-noise ratio greater than 5 (35830 quasars). 

These subsamples are selected according to several different 
criteria. The DR7Q catalogue contains quasars mostly at z < 2 
(see Shen et al. 2011, for their properties) selected with a com¬ 
plex algorithm described in Richards et al. (2002). The DR9Q 
quasar catalogue contains mostly newly discovered quasars, but 
also some DR7 quasars (about 16420) that were re-observed 
during the Baryon Oscillation Spectroscopic Survey (BOSS, 
Schlegel et al. 2007). All quasars have absolute i-band mag¬ 
nitudes brighter than -22, contain at least one broad emission 
line (with FWHM > 500 km s ^) or - if they do not - have in¬ 
teresting/complex absorption features, and have highly reliable 
redshifts. The threshold in absolute magnitude was relaxed for 
DR9Q (M, > -20.5) because BOSS were targeting mostly z > 
2.15 quasars. 

We then searched in the Herschel archive for PACS and/or 
SPIRE images covering a portion of a circle of five arcmin of 
radius around the quasars position for all the ~79000 quasars. 
We searched for the nearest source within a circular region of 
20" centred on SDSS quasar position. We found 224 quasars 
observed by Herschel in point source photometry mode provid¬ 
ing the best instrumental sensitivity and resolution available (114 
from DR7Q, 92 from DR9Q, and 18 from DRIO data archive). 
In Fig. 1 we show the distribution of the redshift values for the 
quasars in our sample. 

3. Data analysis 

3.1. Herschel photometry 

The Herschel Space Observatory is equipped with a 3.5 m tele¬ 
scope that has performed photometry and spectroscopy in ap¬ 
proximately the 55-671 pm range, mapping the “warm” dust 
heated by young stars. There are two photometric instruments, 
the Photodetecting Array Camera and Spectrometer (PACS), 
which observes in three colours (blue 60-85 pm, green 85-130 
pm, and red 130-210 pm), and the Spectral and Photometric 
Imaging Receiver (SPIRE), which observes at 250, 350, and 500 
yum (Poglitsch et al. 2010, Grifhn et al. 2010). 

We perform aperture photometry on PACS and SPIRE 
pipeline-processed images using the annularSkySperturePho- 
tometry task in HIPE version 12.0. We adopt a circular tar¬ 
get aperture centred on the position of the target of radius 10" 
and 15" for blue/green and red PACS images, respectively, and 
22", 30", and 42" for the three SPIRE bands (see Fig. 2). The 
background annulus is chosen between 20" and 40" for PACS 
blue/green bands, between 25" and 45" for PACS red band, and 
between 65" and 95" for SPIRE images (see PACS and SPIRE 
data reduction guide' for details on the PSF shape and the sug¬ 
gested photometry aperture). The default algorithm for the sky 
estimation is the same used in the IDE daophot package, which 
eliminate the outlier pixels, i.e. values with a low probability 
given a Gaussian with specified average and sigma. We define 
the luminosity errors cr^rc as the square root of the quadratic sum 

' http://herschel.esac.esa.int/hcss-doc-14.0/ 


of two terms, ^jNsrc_pixO'sky, and the error on the mean of sky 
jNsrc nix (however the second term is irrelevant because 

' yN,i, 

the sky is measured on a much larger area with respect to the 
extraction region); cTsrc is therefore given by the uncertainty of 
the background times the square root of the number of pixels in 
the extraction region, Nsrc_pix- 

If we do not reach a signal-to-noise ratio (S/N, defined as the 
ratio of the flux density S over the uncertainty error (Tsrc) of at 
least a factor 3, we estimate an upper limit for the flux density 
at 2cr confidence level as IxcTsrc x -\/Nsrc_pix (in this case CTsrc 
is the standard deviation of the flux density mean values in the 
extraction region). Finally we applied aperture correction using 
the task photApertureCorrectionPointSource. 

The sensitivity of the Herschel instruments depends on many 
factors, but to a first order in all the observing modes it scales 
with the inverse of the square root of the on-source observation 
time. The typical sensitivity for a 5cr detection in one hour of 
observation is of the order of 5 mJy for PACS in the 70/100 
/rmband and 10 mJy for SPIRE in the 160 yumband^. In Fig. 
3, we compare the PACS AND SPIRE fluxes in all the bands 
with this approximate sensitivity threshold. The median expo¬ 
sure time for our data is comparable, and all our upper limits 
are located below the threshold for detectability in one hour of 
observation and the detections (with very few exceptions) are 
located above these lines. 

3.2. Ancillary data 

We build the spectral energy distribution (SED) of the 224 
quasars by adding photometric points from ~ 0.3 to 22 pm band 
to the Herschel data, taking advantage of optical and infrared 
all-sky surveys (SDSS, 2MASS, WISE). 

In the optical range, SDSS photometric system comprises 
five colour bands (u, g , r, i, z) from about 0.3 to 0.9 yum . We 
choose the quantity psfMag, for which the total flux is deter¬ 
mined by fitting a PSF model to the object. 

The Two Micron Sky Survey mapped the entire sky in the 
J (1.24 yum), H (1.66 yum), and K (2.16 yum) near-infrared bands 
with a pixel size of 2.0". We considered sources in 2MASS from 
the All Sky Point Source Catalogue and the Reject Table within 
2" matching radius. If the QSO is not detected, we downloaded 
the image of the field of view around the quasars position from 
the 2MASS archive, and we used the IDE routine aper (adapted 
from DAOPHOT) to estimate the limiting magnitude at the po¬ 
sition of the quasar. In detail we derive the magnitude from a 4" 
radius region centred on the expected source position and we es¬ 
timate the sky background in an annular region with inner radius 
of 14" and outer radius of 20". To decide whether the source is 
detected in the image, we repeat the measurement in 1000 cir¬ 
cular regions randomly selected around 100" from the source. 
We then build the distribution of all 1000 flux values (mostly 
from empty regions of the sky) and if the flux of the source is 
higher than the 95th percentile of this distribution, we consider 
this measurement a detection, otherwise the 95th percentile as a 
2cr upper limit. 

The Wide-held Infrared Survey Explorer (WISE) is an in¬ 
frared space telescope that performed an all-sky astronomical 
survey in the W1 (3.4 pm), W2 (4.6 yum), W3 (12 yum), and W4 
(22 yum) bands. All 224 quasars in our sample are detected in 
at least one WISE band and are present in the WISE All-Sky 
Source Catalog. When the flux measurement has a signal-to- 

^ http://herschel.esac.esa.int/Docs/Herschel/html/ch03s02.html 
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Fig. 3. Left panel: observed source fluxes in the PACS bands at 70/100 /an and 160 /rm (green and red points, respectively), compared with the 
detector sensitivity threshold at 5cr for one hour of observation (shown as dash-dotted green and red lines). Right panel: same figure for the SPIRE 
band at 250 /rm, 350 /rm, 500 /rm (blue, green, and red points, respectively). The three dash-dotted lines are the detector sensitivity thresholds in 
the three bands at 5cr for one-hour observations. 
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Fig. 4. Top panel: Example of SED fitting for three quasars SDSS J015950.24-H002340.8, SDSS J022430.60-000038.9, and SDSS J021434.78- 
004242.6 of different SED and spectra quality. Points in yellow represent the luminosity of quasars (in solar luminosity units) from SDSS, 2MASS, 
WISE, PACS, and SPIRE images plotted versus the rest-frame wavelength (from 0.1 to 1000 /rm). We plot in green the starburst template, in 
orange the torus and the hot black-body model, in brown the old star, and in blue the accretion disk component. The results of the composition is 
the best fit, shown in red. The hatched area represents the uncertainties for the normalization of each model components obtained from the Monte 
Carlo analysis. Bottom panel: example of the fit of the SDSS optical spectra for the same quasars presented above. In blue we plot the different 
components of the emission lines (narrow and broad) and in cyan the Fe emission lines template. In the bottom panel we show the residuals of the 
fit. 


noise ratio less than 2, the magnitudes are replaced with the 2cr 
brightness upper limit in magnitude units. 


To characterize the radio properties of our quasars, we follow 
the same approach adopted by Shen et al. (2011) to include the 
radio loudness parameter in the SDSS-DR7 catalogue quasars. 
The standard radio loudness parameter is defined as the ratio be¬ 
tween the flux density (fy) at rest-frame 6 cm and 2500 A (R = 


fecm/fisooA’ Ji™g al. 2007)^. About 25% of the quasars 
(i.e., 50/224) have been detected in the FIRST catalogue and 
34/224 have a radio loudness parameter greater than 10 (~ 15% 
of our quasars are radio loud). 

^ We match our sample with the VLA FIRST Survey at 20 cm with a 
searching radius of 30" and we derive the rest-frame 6 cm flux density 
from the FIRST integrated flux density at 20 cm assuming a power-law 
slope of tty = -0.5. The rest-frame 2500 A flux density is determined 
from the fit to the continuum around 5100 A fitted with a power-law 
shape (fi = Ad" ). 
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3.3. SED fitting 

The infrared emission is dominated by the black-body emission 
of dust at different temperatures, heated by the optical and UV 
emission of the AGN and of young and old stars. A major diffi¬ 
culty in deriving the star formation rate from infrared luminosity 
is to disentangle contributions from different sources of heating. 
Typically the dusty torus can reach dust sublimation tempera¬ 
tures of the order of 1700 K and the peak of the emission is at 
shorter wavelength (less than 100 gm) than that produced by 
warm dust around young stars. In star-forming galaxies there is 
a considerable overlap around 3-20 gm where the polycyclic aro¬ 
matic hydrocarbon (PAH) emission may dominate the IR emis¬ 
sion. However, this is not problematic because PAH features are 
almost absent in AGNs since molecules should be destroyed by 
the extreme optical-UV and soft X-ray radiation in AGN (Roche 
et al. 1991; Siebenmorgen et al. 2004). At shorter wavelengths 
(around 3-4 gm), the emission is dominated by older stars and 
luminosities are thus strongly correlated with stellar mass (Meidt 
et al. 2012). 

In the literature there are many infrared templates for mod¬ 
elling galaxy infrared spectral energy distribution. Some of these 
models are multi-dimensional, in the sense that there is a grid of 
parameters from which one can derive many physical character¬ 
istics of the system (e.g. Silva et al. 1998, Draine & Li 2007). 
Other are one-dimensional, i.e. they rely on a single parameter 
to characterize the shape of the model (e.g. Chary & Elbaz 2001 
and Rieke et al. 2009). 

The measurements in low IR-luminosity quasars is of fun¬ 
damental importance in this analysis since the sources that are 
not detected in Herschel images are likely the best candidates 
to reveal the quenching effect of powerful AGN winds on the 
star formation. Unfortunately, there are no public SED fitting 
codes that can properly deal with upper limit photometric points. 
Therefore we developed our SED fitting code, in which the best 
fit is the result of the sum of different components that describe 
the accretion disk emission, the unabsorbed stellar population, 
the emission of dust heated by the AGNs and by blue young 
stars. Each component is the result of the combination of many 
templates that we specify below. The weight Wj of each template 
can be found by minimizing the function 


a" = s,- 


(Fi - 

AFi 


where E, and AE, are the fluxes and error measured in each pho¬ 
tometric band and E^ is the flux in band i for template j. Since 
this is a linear bounded problem, we determine the weight of 
each template using the idl routine BVLS.pro (bounded vari¬ 
able least-squares) by Michele Cappellari"^. This routine solves 
the linear least-squares problem by finding the weights of each 
template in order to minimize the x^ function. To estimate er¬ 
rors on fitting parameters, we adopt a Monte Carlo procedure. 
Eor each object we generate 1000 realizations of the observed 
SED with random extraction of photometric points normally dis¬ 
tributed around the observed flux with sigma given by the errors. 
The photometric upper limits are treated as null detections with a 
positive error bar at 2cr and we extract random values uniformly 
distributed in this range. Eor each realization of the observed 
SED we computed all the relevant physical quantities, for exam¬ 
ple the luminosity of each component (Lagn^ Lsb, etc.). Einally, 
we estimate the errors as the percentile at 10% and 90% levels 
of the distribution of all the luminosities. If the best model does 

http://www-astro.physics.ox.ac.uk/ mxc/software/ 


not require a starburst component, or the distribution is consis¬ 
tent with a null median value, we consider the 90th percentile of 
the distribution as a 2cr upper limit. 

Eor the SED fitting for each component we adopt the follow¬ 
ing templates: 

- the infrared emission of the starburst component was mod¬ 
elled with a linear combination of the Chary & Elbaz (2001) 
templates. These models provide synthesized SEDs of 105 
average SED templates of local luminous and ultra-luminous 
purely star-forming infrared galaxies (LIRGs and ULIRGs) 
in a sequence of increasing infrared luminosity^; 

- the AGN torus emission was modelled with a linear combi¬ 
nation of DUSTY models by Nenkova et al. (2008). To se¬ 
lect a few likely independent models (from more than 10^), 
we follow Roseboom et al. (2013) who demonstrate that it 
is possible to adequately fit their SED of luminous type 1 
quasars from the SDSS survey with only three models. These 
three models are for tori with inclinations i=0 and i=20 (see 
table 2 of their paper for details); to these we add the models 
at inclinations of 30, 35, 45, 50, 60, 70, and 75 degrees; 

- the emission from three pure black bodies was modelled at 
temperatures 1100, 1300, and 1500 K in order to account 
for the hot dust emission in the near-infrared not adequately 
reproduced by torus models (e.g. Leipski et al. 2014); 

- the stellar emission from an unobscured population was 
modelled with stellar templates from Bruzual & Chariot 
(2003) computed with a Chabrier (2003) IME and metallic- 
ities Z = 0.008,0.02,0.05; we considered models with ages 
5.1 X 10^2.0 X 10^5.0 X lO^l.Ox lO^^, 1.3 x 10'° years; 

- the accretion disk emission was modelled with a linear com¬ 
bination of the four model spectra from figure 1 of Slone & 
Netzer (2012) which are representative of a physical condi¬ 
tions of quasar. 

- for the radio loud quasars it is possible that part of the in¬ 
frared emission is dominated by synchrotron processes orig¬ 
inating in the radio jets. However, the power-law slope of 
the sub-mm/far-infrared continuum is difficult to constrain 
in the infrared band. Eor this reason we prefer not to add a 
power-law component to the fitting templates and highlight 
the radio loud sources in all the plots, remembering that in 
some of these quasars the star formation rate may be overes¬ 
timated. 


In Eig. 4, top panel, we present examples of SED fitting. 

We estimate the star formation rate from the starburst model 
in the rest frame 8-1000 fim adopting a Chabrier initial mass 
function (Chabrier 2003). In detail we use the relation 


SER 

Moyr-i 


= 0.94 X 3.88 X 10-^ 

\ergs ^ 


from Murphy et al. (2011). The factor 0.94 accounts for a slight 
SER underestimation derived using the Chabrier (2003) IME 
with respect to the Kroupa (2001) IME (see Bolzonella et al. 
2010; Pozzetti et al. 2010; Hainline et al. 2011). The stellar 
masses are derived from the normalization of the Bruzual & 
Chariot (2003) templates, since these models are normalized to 
a total mass of one Mq in stars. In the Appendix we test the ac¬ 
curacy of our SED fitting code, in particular with respect to the 
problem of model degeneracy. 

^ We have verified that the choice of the templates does not affect our 
results: adopting the starburst templates by Rieke et al. (2009) or Dale 
et al. (2014), the resultant SB and AGN luminosity differs by less than 
a factor of 3. 
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Fig. 5. Logarithmic flux (in erg s“'cm“^) versus the signal-to-noise 
ratio of the [O III] emission line. We adopt the dotted line as the thresh¬ 
old in flux and S/N to select only spectra with a visible [O III] emission 
line and reliable fit (green points). Red points represent the quasars 
for which we consider the kinematic parameters based on unreliable 
[O III] lines and we do not consider these QSO in the following analy¬ 
sis. 
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Fig. 6. Histogram of the star formation rate distribution 

for all the quasars (top panel) and only for quasars with reliable 
[O III] measurements. The hatched area represent quasars with upper 
limit values on SFRs. 

3.4. Fit of the optical spectra 

In this section we describe the analysis of the spectra used 
to reproduce the profile of the [O III]/14959,5007 and of the 
[O II]/13727,3729 emission line doublet in order to character¬ 
ize the outflow properties. The reduced spectra data used in 
this work are available through the SDSS Data Archive Server 
10 (DAS). The wavelength coverage is Ad ~3800-9200 A with 


a spectral resolution R~1800 - 2200 (Ad ~3600-10500 A and 
R~ 1400-2600 for BOSS spectra). 

We employed a minimization procedure using the IDL 
routine mpflt, Markwardt (2009). Since the optical emission 
lines show complex profiles diverging from a simple Gaussian 
profile, we fit the principal emission lines (i.e. the [O III] dou¬ 
blet, the narrow H/land eventually the He II d4686 with two 
Gaussians, one for a narrow component and one for the wings. 
The intensity ratio is fixed at [OIII] d4959/[0 III]d5007 = 0.33 
(e.g. Dimitrijevic et al. 2007) and the width of the [O III]d4959 
line was fixed to be the same as the [O III]d5007 line and of 
the narrow H^Sline. The velocity offset of the narrow HySline are 
tied to those of the core [O III]d4959,5007 components. The 
[OII] emission line doublet in general is unresolved in our spec¬ 
tra. However, we fit the profile with two emission lines with the 
same cr (each one made of two Gaussian components) with a 
flux ratio of the two components that vary between 1.3 and 1.5. 

The continuum is fitted with a power-law and the broad 
HjS emission line with a broken power-law profile (Nagao et al. 
2006). A major difficulty in the spectral analysis is that the [OIII] 
doublet overlaps with the broad Hy6 emission and optical Fe II 
emission lines (see Fig.4, bottom panel). The Fe II emission 
was modelled with a linear combination of templates, which was 
then convolved with a broadening Gaussian function. The free 
parameters characterizing the Fe II emission are the weights of 
the template, the central velocity, and the velocity dispersion of 
the Gaussian function. The adopted Fe II templates were the 
observed IZwl spectrum (Veron-Cetty et al. 2004). 

In many cases the [O III] and [OII] lines are too weak and 
confused with the continuum and the other fitting components 
to extract reliable kinematic information. Therefore, we adopt 
thresholds in flux and in S/N to discard useless spectra. We only 
consider spectra with 

- S/N>5 for the emission line flux; 

- log(F[o III]) > -11.5 erg s^^ cm^^ 

For the [Oil] line we lower the flux threshold to -12.5 erg 
cm“^ s ', because this line is weaker and not blended with other 
strong features. We exclude an additional 15 sources with high 
S/N and high F[oiii], but with the [OIII]line blended with Fe 
lines or and with a flux that is negligible with respect to the 
other components. The presence of optical Fe emission lines in 
SDSS quasar spectra has been systematically investigated by Hu 
et al. (2008). They found that Fe II emission lines are typically 
redshifted with respect to the [O III] peak, and they suggest that 
the optical Fe emission lines trace gas in the broad line region 
dominated by infall. In Fig. 5 we represent our selection criteria 
on [O III] plotting the flux versus the S/N of the line. Red points 
represent the rejected quasars and green points quasars with bona 
fide clear lines and reliable fits to the profile (confirming by eye 
the quality of our selection criteria). Adopting this approach, we 
define a high quality subsample 1, composed of 125 QSOs with 
trustworthy fitted profiles in [O III] (out of 224 QSOs), and a 
high quality subsample 2, made of 108 QSOs with good qual¬ 
ity fits both in [OIII] and [O II]. These quasars are used in the 
plots involving [OII] in addition to [O III] (we have respectively 
19% and 18% of radio loud quasars in the two subsamples). In 
the following analysis we extract the outflow properties (on a 
case-by-case basis using the [O III], the [OII], or sometimes 
both the line profiles) considering only the quasars that satisfy 
our quality selection criterium. In Fig.6 we plot the histogram 
of the star formation rate in the complete sample and in our high 
quality spectra sample, highlighting the objects with star forma¬ 
tion rate upper limits. The fraction and the distribution of SFR 
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Fig. 7. Left panel: Example of [OIII]/i5007 line decomposition in a narrow and broad component and of different velocity definition used to 
characterize the properties of the outflow. In particular we show velocities at different percentiles of the flux contained in the overall emission line 
profile (from left to right: 2nd; 5th, 50th, 95th, and 98th). Right panel: two times the blue tail, defined as (v50-v02) versus maximum velocity of 
the [O III] line (v98-v02). The dashed line is the bisectrix of the plane and represents the condition for symmetry. Almost all of the objects show a 
broad blue tail. Green points are for radio loud quasars. 
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Fig. 8. Comparison of the [OIII] and [Oil] principal emission line parameters. Left panel: Sigma of the [O III] line versus the sigma of the 
[O II]. The [OIII] lines show broader widths. Right panel: cr excess (defined using the square root of the quadratic difference of the [OIII] and 
[OII] sigma, see text) versus the relative shift of the [O III] peak respect to [OII] peak (both referring to their rest frame position). 


undetected objects is of great interest, since in these objects the 
SF-quenching might be ongoing. In the complete sample 38% of 
the quasars have a star formation rate upper limit measurement. 
This fraction decreases to 26% in the subsample of quasars with 
reliable [O IIII] measurements. We note that although the sensi¬ 
bility of our observations is almost heterogeneous, the distribu¬ 
tion of the IR luminosity for detected and undetected quasars is 
similar and the IR luminosity upper limits are not located in the 
tail of the distribution of detected quasars. 

Owing to the uncertainties in the decomposition of the line 
into a narrow and a broad component, we follow the approach 
of other authors (e.g. Harrison et al. 2014, Perna et al. 2015) 
and we adopt a non-parametric definition for the velocity char¬ 
acterizing the outflow. We reconstructed a synthetic line profile 


adding two Gaussian lines and we define different percentiles of 
the flux contained in the overall emission profile. In particular 
we measure the velocity at the 2nd, 5th, 50th, 95th, and 98th 
percentile, containing respectively 2, 5, 50, 95, and 98% of the 
overall line flux (y02, y05, y50, etc.). The standard deviation (cr) 
of the line profile is the square root of the second moment about 
the mean. In Fig. 7, left panel, we illustrate these definitions 
using an example of a [O III]/15007 line profile. 

In the following analysis, we select the following three pa¬ 
rameters to describe the outflow properties and to explore the 
connection with the SFR; 

- Vbiue, the maximum velocity of the blue wing, defined as 
|vio - V 50 I, linked in some complex way to the energetic of 
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the outflow. In principle, one would infer the mass outflow 
rate, but this measurement is highly uncertain and requires 
the knowledge of the location, mass density, geometry, and 
velocity of the gas outflowing. Assuming simple geomet¬ 
ric arguments, the mass outflow rate is proportional to the 
maximum velocity of the outflowing wind (~ 3vM/R, e.g. 
Maiolino et al. 2012). 

- Avoffset, the velocity offset of the broad wings defined as 
Avoffset = |v'o 5 + ^ 95 l/ 2 -V 5 o. For two well-separated Gaussian 
components, Avoffset is the velocity offset of the broad com¬ 
ponent with respect to V 50 (Harrison et al. 2014). As usual, 
in the literature we interpret the broad, blueshifted emission 
line with respect to the velocity of the [O III] narrow com¬ 
ponent as the result of an outflowing gas moving at larger 
velocity with respect to the quieter large-scale gas. 

- Av„-[o iii]=(cr[o III] ^-<T[o II] ^)°'^, the sigma excess, defined as 
the square root of the quadratic difference between the stan¬ 
dard deviation of the [OIII] and [Oil] lines. This parameter 
estimates the different level of the perturbation of the gas in 
the inner and outer portion of the NLR since, as we discuss 
in Section 4.1, the bulk of the [O III] emission should be pro¬ 
duced in a more inner region than [OII]. 

4. Results 

4 .1 . Comparison of the [O III] and [OII] line profile 

In the literature a blueshifted line seen in absorption is consid¬ 
ered a clear signature of outflow. Instead, the interpretation of 
blueshifted wings seen in emission is ambiguous because differ¬ 
ent scenarios can account for the asymmetry of the line: out¬ 
flows, inflows, or gas in a rotating disk partially obscured by an 
asymmetric distribution of dust. Moreover, if the extinction is 
within the emitting clouds, then outflowing gas would produce a 
line profile with a red wing (Whittle 1985; De Robertis & Shaw 
1990; Veilleux 1991). However, in recent papers a blueshifted 
wing is considered a clear signature of ionized gas moving to¬ 
wards our line of sight (e.g. Bian et al. 2005, Komossa et al. 
2008, Brusa et al. 2015, Villar-Martm et al. 2011a). Here we 
follow the latter interpretation, aware that all our results depends 
on this assumption. 

In Fig.7, right panel, we emphasize the asymmetry of the 
[O III] emission lines plotting Vmax^vgg-voi versus two times the 
blue tail (defined as vso-voi)- The deviation from the bisectrix 
of the plane represents the amount of asymmetry with respect to 
V50. About 80% of the quasars show asymmetric [O III] emission 
lines with significant blue wings. The median deviation from 
symmetry is 211 kms^'but it increases up to 1033 kms^'at 
higher v^n^x- 

In various models of line formation in the NLR, the larger 
widths of the ionized lines is explained by their formation in 
clouds that have higher density or higher orbital or radial ve¬ 
locity, likely because they are closest to the central ionization 
region (De Robertis cfe Osterbrock 1984)®. The larger width of 
the [O III] with respect to [OII] is consistent with the fact that 
lines of higher ionization potential and/or higher critical densi¬ 
ties are produced in inner regions closer to the nucleus. If the 
[Oil] emitting gas likely traces the unperturbed ISM, the ex¬ 
cess of the standard deviation of the [O III] line with respect to 

* The critical density and the ionization potential for [OIII] is 7 10^ 
cm“^ andxev =35.1, higher than the critical density for the [Oil] doublet 
(1.4-3.3 10^ cm“^ eV) and the ionization potential for [Oil] (y„ =13.6 
eV), Nagao et al. (2001). 


[OII] is an estimate of the amount of perturbation of the gas in 
inner regions. 

In Fig.8, left panel, we plot the sigma values derived from 
the analysis of the [OIII] and [Oil] emission line profile. The 
o'jo III] values span a range in velocity between 165 and 950 
kms^*, with a mean value of 405 + 165kms“*. For compar¬ 
ison, Mullaney et al. (2013) in a sample of type 1 quasars found 
a value of the order of 195 km s *(we derived this value combin¬ 
ing their average two-Gaussian component with FWHM of 335 
and 851 kms“^). The value of cr[o iii] > except for a few cases, is 
always larger than ctjq u] . For the [Oil] emission line we find 
o'loii]=185 + 70 kms“\ a value comparable with the stellar 
velocity dispersion, which in quasars is of the order of 200-300 
km s '(depending on the black hole masses). 

In order to characterize the gas outflow velocity with respect 
to the unperturbed gas at galactic scales it is essential to deter¬ 
mine an accurate velocity reference. In our quasar sample the 
stellar absorption-line features cannot be detected, because the 
continuum is dominated by the AGN. The velocity of the core 
of the [O III] line has been used to estimate the systemic redshift 
of the galaxy (Villar-Martm et al. 2011b) and the width of the 
line has been considered a reasonable substitute for the stellar 
velocity dispersion of the galaxy bulge (after properly removing 
the asymmetric blue wings of the line, Greene & Ho 2005). As¬ 
suming that the outer portion of the NLR gas follows pure grav¬ 
itational motions, we test whether the velocity of the [Oil] line 
can be used as a reference velocity for large scale gas. 

In Fig.8, right panel, on the x-axis we show the shift in veloc¬ 
ity of the [O III] peak flux density with respect to the [OII] peak 
both referring to the rest frame wavelength . We measure a mean 
shift of the two peaks of < VsUft >=-20+100km s^'. Unfortu¬ 
nately, the [Oil] doublet is unresolved in the SDSS spectra. The 
ratio of the two lines depends on the electron density of the gas 
and it is between 0.35 and 1.5 (Pradhan et al. 2006). This in¬ 
troduces an uncertainty in our [Oil] velocity reference of 2A, 
corresponding to ~ 80 km s '. For most of the quasars we find 
relative shifts of the two peaks in this range. However, a few 
quasars (~ 20 % of the total) show velocity shifts larger than 80 
kms”'. For comparison, on the y-axis we show Avo-[o mj. The 
quasars that show larger [O III] shift with respect to [OII] also 
have larger Avo-[o iii] values. We believe that this approach is 
promising; unfortunately, the [OII] doublet is unresolved in our 
spectra and so we prefer not to use the peak of this line as a ref¬ 
erence frame and we choose other parameters to characterize the 
outflow properties. 

4.2. SB and AGN dominated quasars 

The total infrared luminosity between 8-1000 pm is the sum 
of different contributions: the torus and dust heated by newly 
formed stars. We divide our sample into two groups according 
to the dominant processes that produce the infrared-luminosity. 
The first group is composed of starburst dominated objects for 
which the starburst-related luminosity dominates the total in¬ 
frared (8-1000 pm) luminosity (LsB/Ttot >0.5) and the sec¬ 
ond group is composed of AGN component dominated objects 
(LsB/Ttot <0.5) (or quasars with upper limits in SFR). In this 
section we explore the properties of these two groups (hereafter 
SFd-QSO and AGNd-QSO). 

We measured the star formation rate from SED fitting in the 
infrared band; however, the [O II] emission line, corrected for 
dust extinction, is another widely used estimator of the SFR. 
Deriving the SFR from the [OII] luminosity is quite a challeng- 
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Fig. 9. Upper panel: Logarithm of the infrared luminosity in the 
8-1000 yum range for the starburst component derived from the SED 
fitting versus the logarithm of the [Oil] luminosity, another estimator 
for the SFR. Bottom panel: Logarithm of the infrared luminosity in 
the 8-1000 yum range for the hot dust component derived from the SED 
model (likely produced by a dusty torus warm up by the AGN) versus 
the logarithm of the AGN bolometric luminosity. Blue and red points 
represent quasars dominated in the infrared by star formation or by the 
AGN emission, respectively. 


ing task because in type 1 quasars the [OII] luminosity is con¬ 
taminated by the emission produced in the NLR. Moreover, the 
conversion from luminosity to SFR has a complex dependence 
on metallicity (Kewley et al. 2004). In Fig.9, upper panel, we 
compare the [Oil] versus the infrared luminosity in the 8-1000 
/tm derived from SED htting. The [OII] luminosity is not cor¬ 
rected for dust absorption because only about half of our spectra 
(z ^ 0.4) cover the Balmer decrement (Ha/Hy6) typically used 
to make the correction. As expected, the dispersion is larger if 
we consider all the quasars, but it is significantly reduced if we 
consider only the SFd-quasars. The red points (AGNd-quasars, 
102 objects) may have L[oii] signihcantly contaminated by the 
AGN, and in fact they show an excess of [Oil] luminosity with 
respect to the blue points (SFd-dominated quasars, 46 objects). 


47 



41 42 43 44 

log L[OIII] [erg s ’] 


Fig. 10. Logarithm of the star formation rate versus the logarithm of 
the [OIII] luminosity. The dashed line shows the correlation by Netzer 
(2009) of slope 0.8. The normalization is from Shao et al. (2010) divid¬ 
ing the AGN bolometric factor by a factor of 3500. The histograms at 
the right and bottom represent the distribution of [OIII] luminosity and 
of star formation rates. 

In the bottom panel of Fig.9 we compare the 8-1000 /urn lu¬ 
minosity for the AGN component versus the AGN bolometric 
luminosity (derived from L[o mj as Lbol = 3500xL[o mj, Heck¬ 
man et al. 2004). The ratio between these two quantities pro¬ 
vides an estimate of the torus covering factor. We And a mean 
logarithmic value of <-l+ 0.5, consistent with the estimate of - 
0.4+ 0.2 of Roseboom et al. (2013) measured in a sample of type 
1 quasars (our [O III] luminosity measurements are not corrected 
for dust absorption and therefore we consider them as lower limit 
values). In the AGNd-quasars sample (86 objects), the L[o mj 
correlates with the infrared emission derived from torus models 
(we And a linear correlation coefAcient of r=0.6). Instead we do 
not And a relation in the SFd-quasars sample (39 objects, r=0.2). 

In Fig. 10 we investigate the relation between the AGN 
[O III] luminosity and the star formation rate of the host galaxy 
of the two groups. There is a well-known close correlation be¬ 
tween the AGN bolomehic luminosity and the star formation 
rate of galaxies (e.g. Netzer (2009)). This relation is usually 
interpreted as a probe of the coupling between AGN growth and 
the star formation by an evolutionary mechanism, e.g. the merg¬ 
ing. Our quasars recover the correlation found by Netzer (2009) 
Lacn ^-'SFR that links the AGN luminosity to the SFR. 

4.3. SB and outflow properties 

We now consider possible relations between the SFR and the 
kinematic parameter characterizing the outflowing gas and we 
discuss these results in a general context in Section 6. We fo¬ 
cus flrst on the SFd-QSO sample (Fig.11, top panel). We And 
an upward trend between the logarithm of the SFR and the three 
kinematic parameters considered, but no signiflcant correlation. 
For the AGNd-QSO, we do not And a trend between the outflow 
properties and the SFR. In Fig. 11, bottom panel, we consider the 
SFR versus the three same kinematic parameters considered in 
the top panel. Even separating the radio loud and radio quiet 
quasars we did not And a relation between the SER and kinemat¬ 
ics of the ionized gas, as expected if the dominant source pow- 
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Weak-Out. QSO 

Strong-Out. QSO 

zbin 

Det/ul 

<SFR> 

Det/ul 

<SFR> 

^blue 

0 -0.2 

12/1 

1.23(1.07) 

12/0 

1.41 (1.44) 

0.2-0.4 

16/7 

1.27(1.29) 

13/3 

1.69(1.51) 

0.4-0.6 

9/6 

1.32(1.37) 

8/4 

1.44(1.31) 

0 .6-1 

9/7 

1.89(2.01) 

13/5 

2.25 (2.15) 

^^ofFset 

0 -0.2 

12/1 

1.25(1.10) 

12/0 

1.39(1.41) 

0.2-0.4 

16/6 

1.27(1.29) 

13/4 

1.67(1.50) 

0.4-0.6 

11/5 

1.34(1.40) 

6/5 

1.43(1.26) 

0 .6-1 

7/6 

1.95(1.97) 

15/6 

2.20(2.15) 


Avo-[o III] 


0 -0.2 

10/1 

1.22(1.07) 

12/0 

1.45(1.50) 

0.2-0.4 

13/5 

1.42(1.34) 

11/4 

1.63(1.55) 

0.4-0.6 

9/5 

1.44(1.36) 

6/5 

1.29(1.28) 

0 .6-1 

7/5 

1.93(1.95) 

10/5 

2.18(2.09) 


Table 1. Statistic results for the SFR plotted in bin of redshift (see Fig. 
13). For the three kinematic parameters adopted, we report in each bin 
of redshift (col. 1) the number of the detection and upper limits and 
the mean SFR in logarithmic scale (in parenthesis the value obtained 
excluding radio loud quasars) for unperturbed quasars (col. 2 and 3) 
and outflow dominated quasars (col. 4 and 5). 


ering the outflow are supernovae winds or the AGN. We then 
searched for a connection between the kinematics of the ionized 
gas and the properties of the AGN (e.g. L[oiii], itot(8-iooo/jm), 
^Edd), but again we do not find any linear relations. 

4.4. Testing the negative AGN feedback modei 

The SFR of star-forming galaxies correlates with the stellar mass 
in the local Universe, as well as in the distant Universe out to z ~ 
2 (e.g. Noeske et al. 2007; Daddi et al. 2007; Elbaz et al. 2007). 
This close correlation (dispersion of ~ 0.3 dex) is called the main 
sequence and it can be used to define the outliers, i.e. quenched 
and starburst galaxies. Later on, the specific star formation rate 
(sSFR), defined as the SFR per unit stellar mass, is used to in¬ 
vestigate the dependence of the SFR with many parameters char¬ 
acterizing the galaxies or the environment. For example, a de¬ 
pendence of the sSFR on stellar mass and redshift (Bauer et al. 
2005) was discovered. Galaxies at redshift equal to one tend 
to have star formation rates that are a factor of ten higher than 
the local ones. Also, high redshift galaxies are typically more 
massive than low-redshift star-forming galaxies (Seymour et al. 
2008). 

In this work, we compute the SFR in redshift bins inde¬ 
pendently for strong and weak outflow galaxies (assuming sim¬ 
ilar stellar masses in each redshift bin for the two groups). 
In order to divide the quasars into the two groups, we con¬ 
sider the distribution of Vbiue, Avoffset and Avo-[o mj and as the 
threshold we choose the median of the respective distributions 
(Vbiue>556, Avoffset >166 and Av^p iii]>335 kms^'), see Fig. 12. 
In Fig. 13 we show the logarithm of the SFR versus redshift plot¬ 
ting strong and weak outflow quasars according to these three 
criteria. We have separated the quasars into four redshift bins 
(z<0.2, 0.2<z<0.4, 0.4<z<0.6, and z> 0.6, see Table 1) and in 
each redshift bin we compare the mean star formation rate for 
the two classes separately, taking error bars and upper limits in 
the SFR measurements into account. 

It is not an easy task to calculate a weighted mean with up¬ 
per limits in a dataset. We adopt the following strategy. For each 


Variables 

Disp. 

r(N) 

Sign. 1 

Agn dominated quasars 

log SFR vs Max velocity 

0.2 

0.44(39) 

95% 

SFR vs Wings vel. offset 

0.2 

0.40(39) 

95% 

log SFR vs (o-[OIII] - o-[OII]) 

0.2 

0.45(33) 

95% 

SB dominated quasars 

log SFR vs Max velocity 

0.5 

0.23(86) 

NO 

log SFR vs Wings offset 

0.5 

0 .11(86) 

NO 

log SFR vs A VSIG[OIII] 

0.5 

0.22(75) 

NO 


Table 2. Col. 2: The variance of the intrinsic scatter. Col3:the linear 
correlation coefficient between the dependent and independent variables 
(Pearson product-moment correlation coefficient). The maximum pos¬ 
itive (negative) correlation is indicated by a correlation coefficient of 1 
(-1). An absolute value between 0.4 and 0.7 could indicate moderate 
association. Col.4: Significance of the linear correlation. A value of P 
less than 0.01 (0.05) imply a significant measurement at 1%(5%), id.e. 
there is 1% (5%) percentage chance that the finding is non-significant. 
N is the size of the sample. 


quasar we consider the histogram distribution of SB luminosity 
(i.e. the distribution of the luminosity measurements obtained 
through a Monte Carlo analysis). Some of the QSO are consid¬ 
ered upper limits and their distribution peaks at zero luminos¬ 
ity, others are detections and their distribution resemble a Gaus¬ 
sian distribution. We randomly extract lO'* SB luminosities from 
the distributions of all quasars separating them in four classes of 
redshift and in strong and weak outflow quasars. In these four 
redshift bins, we obtain two histograms of the total luminosities 
that take into account detected and undetected sources. We eval¬ 
uate the mean and the error at the 10% and 90% percentiles of 
the distribution and we overplot the median points in the figure. 
Since we may have overestimated the star formation rate for ra¬ 
dio loud quasars, and the type of the AGN feedback on the host 
could be different, we repeat the same analysis excluding radio 
loud quasars (Fig. 13, bottom panel). 

We note an increase in the mean star formation rate with red¬ 
shift with a similar slope in the two groups. This is expected be¬ 
cause it is well known that the cosmic star formation rate peaks 
around z~ 2-3, after which it decreases by an order of magnitude 
to the present Universe (e.g. Hopkins et al. 2006). More impor¬ 
tantly, we find that the star formation rate is similar in the two 
groups of quasars and, typically, strong-outflow quasars show 
values only slightly larger than weak-outflow quasars, even ex¬ 
cluding the radio loud quasars. In each redshift bin the discrep¬ 
ancies between the median values in the two groups are slightly 
reduced if radio loud galaxies are excluded from the analysis. 
We have found no evidence that the SF in the host is suppressed 
in the presence of strong outflows comparing the SFR in bin of 
redshift in quasars characterized by strong or weak outflow sig¬ 
natures. 


5. Main sequence of quasar host galaxies. 

One of the best ways to take into account the effects due to dif¬ 
ferent galaxy masses is to normalize the SFR to the stellar mass 
obtaining the so-called specific star formation rate (sSFR). How¬ 
ever, given the difficulties in measuring stellar masses in quasar 
host galaxies from SED fitting in the IR band caused by the 
strong AGN contamination, we adopt this approach as a con¬ 
sistency check. The offset from the SF main sequence provides 
evidence that the SF in the galaxy is enhanced or quenched with 
respect to most of the galaxies. Therefore, we derive the stel- 
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Fig. 11. Relation between the SFR and the three parameters characterizing the outflow (vbiue. Avoifset, Av£r[o iii]) for the starburst dominated 
quasars (blue points, top panel) and the AGN dominated quasars (red points, bottom panel). There is no clear relation between outflow and SF 
considering the AGN and SB dominated quasars. If we consider only the SFd-dominated objects (blue points) a positive trend emerges, although 
with a weak significance and a large scatter. 



Fig. 12. Histogram of the distribution of the three parameters characterizing the outflow. The red dashed line marks the median of the distribution. 
From left to right: vyue. Avojfset, Avo-[o iii]. The green hatched area represents the radio loud quasars. 


lar mass from the normalization of the Bruzal & Chariot stellar 
template. 

As mentioned above, measuring stellar masses in type 1 
AGNs and in quasars in particular is extremely challenging. 
However, in many cases the near IR SEDs of our sources around 
1-3 yum are dominated by the emission of old stars, allowing a 
reliable estimate of the stellar mass. We consider convincing the 


stellar mass derived from SED fitting in which the luminosity 
of the old star components contribute at least 25% to the total 
luminosity in the 1-3/rm range. We show in Eig.l4, right panel, 
the distribution of the stellar mass values obtained for the whole 
sample of quasars (the dashed area represents uncertain mea¬ 
surements, about 30% of the total measurements). By consid¬ 
ering only the reliable values we obtain a mean value of 10'°'^ 
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Fig. 13. Upper panel: Star formation rate plotted in four bins of redshift for the complete sample of quasars. We adopt different criteria for 
separating strong-outflow (small magenta points) and weak-outflow galaxies (small green points), respectively higher and lower than the median 
value of the distrihution of the three parameters (vbiueffirst panel), Avoffset (second panel), and Av£r[o iii] )• Bottom panel: same as top, hut excluding 
radio loud quasars. 




Redshift 


Fig. 14. Left panel: Histogram of the logarithm of the stellar masses in Mq derived for our sample of type 1 quasars. The hatched area represents 
uncertain values. Right panel: Evolution of the main sequence SFR with redshift. The solid line represents the relation between sSFR and redshift 
and the two dashed lines are a factor of two above and below this line. 


Mq, consistent with the range found by Mainieri et al. (2011) 
between 8x10^ and 10 *' Mq for a sample of type 2 quasars from 
the COSMOS survey. 

In Fig. 14, right panel, we plot the sSFR versus redshift 
and we compare these measurements with the main sequence 
(adopting the Elbaz et al. (2011) formulation^). The sSFR of 
the quasars are in agreement with the sS FRms relation found for 
normal star-forming galaxies, although with a very large scatter 


’’ sS FRms = 26 x in [Gyr *] with ^cosmic Gyr 


that is due in part to the uncertainties on the stellar mass mea¬ 
surement. We investigate whether there is a relation between the 
position in this plane and the kinematic properties of the out¬ 
flow. In Fig. 15, we compare the offset from the main sequence 
of SF (sSFRoffset. a measure of the starburstiness of quasars) 
to the three kinematic parameters. Again, analogously to the 
SFR, we do not find a clear relation between these quantities 
and sSFRoffset- 
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Fig. 15. sSFR offset from the main sequence versus the three emission line parameters characterizing the outflow: Vbiue, Avoifset, Avo-[o iii]- The 
green points represent the radio loud quasars. 


6. Discussion 

6.1. Evidence fora SB-outflow connection? 

The origin of the fast ionized gas outflows is controversial since 
in general it is not possible to establish unambiguously whether 
they are produced by the intense star formation, by the AGNs, 
or by a combination of the two processes. Galactic winds pro¬ 
duced by outflows driven by thermal energy from supernova ex¬ 
plosions or by stellar winds typically exhibit outflow velocities in 
the 100-500 kms * range (Rupke et al. 2005) and outflows with 
significantly higher velocities (|v| > 1000 kms^*) are generally 
interpreted as driven by an AGN (e.g. Tremonti et al. 2007, Ci- 
cone et al. 2014). However, according to Diamond-Stanic et al. 
(2012) ram pressure from supernovae and stellar winds is suffi¬ 
cient to produce outflows with velocity higher than 1000 km s * 
without needing to invoke feedback from an active galactic nu¬ 
cleus. 

Blue star-forming galaxies show a positive correlation be¬ 
tween SFR and outflows. For example, Weiner et al. (2009) 
found that the outflow velocity scales linearly with SFR as 
'^wind ~ SFR**-^ in blue star-forming galaxies at z ~ 1.4, similar 
to the scaling found by Martin (2005) in low-redshift ULIRGs. 
For quasars no obvious relation between the wind velocity and 
infrared luminosity has been discovered. For example, Veilleux 
et al. (2013) found that the velocity of the OH absorption line 
does not correlate with the star formation rate in a sample of 
ULIRGs and QSOs. A similar result is found by Rupke & 
Veilleux (2013), considering both ionized and neutral gas out¬ 
flow in six ULIRGs. These authors did not find any rela¬ 
tion between the OH outflow velocities and the star formation. 
However, they And a non-linear relation between the velocity 
of the molecular outflow and the AGN luminosity, e.g. more 
blueshifted OH components are found in more AGN luminous 
sources. 

We have looked for a relationship between the SFR and the 
strength of the outflow, considering many different outflow pa¬ 
rameters (e.g. Vbiue. Avoffset, Avo.[o III])- We consider separately 
the two classes of SFd and AGNd quasars (Fig. 11, bottom 
panel), in which the driver mechanism of the outflow is prob¬ 
ably different. No clear trend emerges between the star forma¬ 
tion rate and the velocity of the outflow. However, we do not 
And any relation between the outflow and the properties of the 
AGNd quasars, as is expected if the outflow is ultimately driven 
by the radiative output of the quasar and as was found instead 


for molecular outflows in quasars by Veilleux et al. (2013). We 
cannot rule out that both AGN and stellar processes contribute 
to driving the outflow and weaken any relations, especially in 
quasars with a high star-forming activity. This is in agreement 
with what has been found by other previously cited authors, i.e. 
we do not And any obvious relation between the wind velocity 
and infrared luminosity in the quasar sample. 

If we separate the quasars into two classes of quasars with 
weak and strong outflow, an interesting result emerges. We do 
not And that the mean star formation rate is lower in host galax¬ 
ies that exhibit the strongest outflows, as predicted by the basic 
AGN negative feedback model; instead, it is consistent or even 
higher than in weak outflow quasars. However, it is worth noting 
some important shortcomings and uncertainties of this approach. 

- There could be a significant delay between the starburst 
phase and the peak of nuclear optical AGN activity (e.g. 
Matsuoka et al. 2011, Davies et al. 2007; Bennert et al. 
2008; Schawinski et al. 2009; Wild et al. 2010, Yesuf et al. 
2014) estimated in the range of 100-400 Myr. The quasar 
phase is expected to last for a typical lifetime of 16-25 Myr 
(Gon 9 alves et al. 2008) consistent with the estimate of > 10^ 
yr by Jakobsen et al. 2003. The timescale for star formation, 
i.e. over which star formation remains constant, can vary by 
a few 100 Myr (e.g. de Grijs 2001). Therefore, a substan¬ 
tial time lag between the quasars phase and the quenching of 
the star formation rate in the host galaxy could wash out the 
most significant signatures of negative feedback. However, 
we note that if powerful winds are able to sweep away the 
interstellar dust, we would expect an underluminous infrared 
galaxy host, even if the star formation proceeds. 

- Another possible interpretation of this result is that the AGN 
is stimulating the SFR as predicted by a positive feedback 
scenario. In this interpretation the AGN outflow triggers star 
formation in the gas-rich galaxies by overcompressing cold 
dense gas and thus providing positive feedback. AGNs with 
pronounced radio jets exhibit a much higher star formation 
rate than the purely X-ray selected ones (Zinn et al. 2013). In 
another case (Cresci et al. 2015), both types of feedback are 
observed in the same galaxy: the outflow removes the gas 
from the host galaxy (negative feedback), but also triggers 
star formation inducing pressure at the edges of the outflow 
(positive feedback). We expect to see the positive feedback 
effect especially in radio loud quasars. Interestingly, exclud- 
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ing radio loud quasars in our analysis, the difference between 
the mean SFR in the two groups decreases. 

- The blueshifted wings of [O III] are typically seen only in the 
inner kpc of NLRs in local AGNs. So their detections imply 
outflows over regions that are much smaller than the star¬ 
forming portion of typical host galaxies and therefore these 
outflows may not be able to quench the star formation rate. 
Integral field spectroscopic data are crucial to answering to 
this question. 

- Since the Herschel spatial resolution does not allow us to 
measure the SFR in the nuclear region (on a subkilo parsec 
scale), we cannot rule out the possibility that the impact of 
negative feedback even in the most luminous systems only 
has an affect in the very central regions (below 100 pc), as 
suggested recently in the model of Roos et al. (2015). 

- In this work we have considered the SFR, since the stellar 
mass in type 1 quasars is a difficult parameter to determine. 
It is possible that the sSFR correlates with the mass loading 
factor or the kinetic power of the outflow, not with the veloc¬ 
ity. However, we cannot measure it with this data so we use 
the velocity as a proxy of the strength of the outflow. 

Keeping in mind the discussed caveats, our results seem 
to disfavour the AGN negative feedback scenario, according to 
which AGN outflows on the QSO phase sweep out the gas from 
the host galaxy suppressing the star formation. 

Our results are more in agreement with works that seems to 
demonstrate that AGNs do not regulate starbursts in the overall 
galaxy. For example Debuhr et al. (2010), found that SMBH 
feedback has little effect on the number of stars formed. In the 
theoretical model of Cen (2012) starburst and AGN growth are 
not coeval in this model and AGN activities are expected to out¬ 
live the starburst, in agreement with observations (e.g. Geor- 
gakakis et al. 2008). In a recent simulation Gabor (& Dave (2014) 
found that AGN feedback drives bursty, high velocity, hot out¬ 
flows with mass rates peaking briefly near the SFR but that have 
little impact on the star-forming gas disk. The authors conclude 
that the star formation is not quenched by AGN feedback, and 
they suggested that the key physical element required to produce 
a galaxy red sequence is not expelling gas from galaxies, but 
rather preventing gas from accreting. Other processes than AGN 
negative feedback are able to transform a star-forming galaxy 
into a quiescent system, for example the so-called strangulation 
mechanism (e.g Larson et al. 1980, Balogh et al. 2000). In this 
case star formation in the host continues until the available gas 
is completely used up and the star formation is not suddenly 
quenched as in the negative feedback model. In a recent analysis 
Peng et al. 2015 found that strangulation is the primary mecha¬ 
nism acting in local galaxies. 

The increased probability of an AGN being hosted by a star¬ 
forming galaxy (e.g. Santini et al. 2012) may be a consequence 
of the relationship between gas content and AGN activity. If this 
is the case, the relationship between AGN bolometric luminos¬ 
ity and the SFR found by Netzer (2009) naturally arises. This 
relationship might ultimately be caused by secular mechanisms 
of gas inflow or merger mechanisms between galaxies. What 
emerges clearly is that the complex interplay between AGN out¬ 
flows and star formation rate is far from being understood. 

7. Summary and conclusions 

The aim of this paper is to test the negative AGN feedback sce¬ 
nario in its simplest version according to which QSO driven out¬ 
flows are able to suddenly quench the SFR in the host, clearing 


out or heating the interstellar medium of the galaxy. We have 
considered 132 quasars from the SDSS quasars observed in pho¬ 
tometric mode by the infrared satellite Herschel at redshifts of 
less than 1 (to cover the [O III] line in the SDSS spectra) with 
high quality fits and spectra. 

To investigate the presence of an outflow, we focus on the 
[O III] and [OII] emission line, in particular on three parameters 
(the maximum blue velocity, Vbiue=kl0-v50|; the offset of the 
broad emission, Avoffset =|v95h-v5|/2-v50; and the cr excess with 
respect to the [Oil] line, Av^p iii]=(o-([OIII] )2 - cr([OII] )2)0-5) 
assuming that a broad, blueshifted wing is a signature of an out¬ 
flow. 

Since the SFR evolves with redshift and it is dependent on 
the stellar mass of the galaxy (galaxies at higher redshift are typ¬ 
ically more massive and more star forming), we consider the 
mean SFR in bins of redshift for strong and weak outflow galax¬ 
ies. We do not see that strong outflow galaxies have a SFR lower 
than weak-outflow galaxies as was predicted by the negative 
AGN feedback scenario. Instead, we found that the SFRs in the 
two groups are comparable or even higher for the outflow dom¬ 
inated galaxies. We found an analogous result considering the 
offset from the main sequence of the sSFR. The stellar masses 
derived from SED fitting in type 1 quasars is uncertain because 
of the AGN contamination in IR. However, for the quasars with 
trustable measurements in which the stellar component domi¬ 
nates in the range 1-3 fim, we find that the redshift evolution of 
the specific SFR for the hosts in our sample of type 1 QSOs is in 
agreement with that observed for star-forming galaxies. 

To investigate the main driver of the outflow, we divided our 
sample into two groups: the AGN-dominated quasars (defined 
as LsB/Ltot <0.5) and the SF-dominated quasars (LsB/Ltot >0.5). 
We find that by focusing only on SF-dominated quasars we did 
not find a correlation between the SFR and the velocity, as found 
for starburst galaxies; rather, a slight positive trend was found. 
Instead in the AGN-dominated quasars do not find any clear rela¬ 
tionship between the outflow velocities and the SFR or the AGN 
properties (Eddington ratio, [O III] luminosity, or total infrared 
luminosity). 

We conclude that the negative feedback scenario is dis¬ 
favoured by our results and that other possible mechanisms could 
be responsible for quenching star formation. In our large sam¬ 
ple of quasars, a high level of SER coexists with powerful out¬ 
flow and, moreover, outflow dominated host galaxies have star 
formation rates consistent with or slightly larger than galaxies 
showing weak outflows. Different possibilities can explain the 
lack of observational evidence of negative feedback, including 
the possibility that AGN driven outflows are not the main mech¬ 
anism by which the star formation in blue star-forming galaxies 
is suppressed. 
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Appendix A: 

We explored the accuracy of our SED fitting code, since the best 
model that reproduces the photometric data points can poten¬ 
tially be subject to model degeneracy. To evaluate the uncertain¬ 
ties on the main quantities derived from the fit, i.e. Lagn, Tsf, 
and Mstar, in relation to the photometric uncertainties and model 
degeneracy, we performed some tests. 

First of all, we built various theoretical SEDs, extracting 
photometric fluxes from our best fitting models. In each case 
the integrated luminosity of the SF and AGN component were 
known. 

Second, we reproduced real cases, adding observed, realistic 
random errors to the theoretical photometric points (substituting 
randomly some photometric points with upper limits). Then we 
fit this SED obtaining Lsf Lagn observed. We repeated this pro¬ 
cedure for 100 times. 

In Fig.Al, left panel, we show an example of a result of this 
analysis for one theoretical SED. We plotted the distribution of 
the values obtained in 100 iterations, and we compared the ob¬ 
tained values with the theoretical value shown with a dashed line. 
In each iteration the best model is composed of different tem¬ 
plates; however, the distribution of the integrated luminosities 
are centred on the expected value. In Fig.Al, right panel, we 
report in percentage the distribution of the relative errors for the 
main quantities Lagn, Lsf, and Mstar in all the iterations com¬ 
pared to the errors that we obtained in the theoretical case. 

As a consistency check, we compare the colour obtained by 
Veilleux et al. (2009) for a sample of ultraluminous infrared 
galaxies and Palomar Green quasars with Spitzer. We derive 
the flux ratios at different wavelengths from the SEDs best fit¬ 
ting model. In Fig.A.2, left panel, we compare L'so^jm/L'is^m 
with respect to the flux ratio between 5 and 25 /rm(mid-IR) 
and between 40 and 120 /im(far-IR). In the right panel, we 
plot with respect to In both cases our 

quasars extend the relation found by Veilleux et al. (2009) down¬ 
ward. 
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Fig. A.l. Example of a SED fitting test. Left panel: Distribution of the AGN luminosity, SB luminosity, and stellar mass obtained in 100 iterations 
obtained from theoretical SEDs, and by varying the photometric points according to the observed uncertainties compared to the theoretical value 
represented by the dashed line. Right panel: as before, but we plot the distribution in percentage of the relative errors for the parameters in each 
iteration. 
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Fig. A.2. Comparison of SED colours derived from the SED with respect to a sample of PG quasars with Spitzer measurements and ULIR 
galaxies (black points, from Veilleux et al. 2009) and our sample of quasars (magenta points). 
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